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ABSTRACT 

It  is  shown  that  when  an  acceleration  wave  propagates  in  a  hyperelastic 
rod  with  slowly  varying  cross-section,  the  transport  equation  for  the  wave 
intensity  is  a  generalized  Riccati  equation.  The  three  coefficients  in  the 
equation  all  depend  on  the  material  properties,  but  only  the  coefficient  of 
the  quadratic  term  is  independent  of  the  effect  of  area  change.  Three 
theorems  are  proved,  based  on  the  use  of  comparison  equations,  which  establish 
that  in  general  the  acceleration  wave  intensity  will  become  infinite  (escape) 
after  the  wave  has  propagated  only  a  finite  distance  along  the  rod.  The 
existence  of  thresholds  for  the  initial  intensity  are  also  established  in 
certain  cases,  with  their  most  notable  property  being  that  as  the  initial 
intensity  decreases  towards  the  threshold,  so  the  distance  the  wave  propagates 
to  escape  increases  without  bound. 
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SIGNIFICANCE  AND  EXPLANATION 


This  work  studies  the  propagation  of  an  acceleration  wave  in  the  class  of 
nonlinear  elastic  materials  known  as  hyperelastic  •  These  are  materials  in 
which  the  stress  potential  depends  only  on  the  displacement  gradient.  In 
particular/  the  effect  of  wave  propagation  along  a  rod  of  slowly  varying 
cross-section  is  studied,  since  in  many  physical  situations  sucii  rods  are 
often  to  be  found.  The  assumption  of  a  slowly  varying  cross-section  allows 
the  problem  to  be  approximated  by  a  one- dimensional  situation,  with  the  area 
variation  either  modifying  terms  or  introducing  new  ones. 

The  equation  governing  the  acceleration  wave  intensity,  the  transport 
equation,  is  shown  not  to  be  the  Bernoulli  equation,  which  usually  arises,  but 
the  generalized  Riccati  equation.  In  the  case  of  a  rod  of  constant  cross- 
section  this  reduces  to  a  degenerate  Bernoulli  equation  which  may  be  solved 
exactly  to  yield  results  which,  in  conjunction  with  experiment,  permit  the 
determination  of  the  material  characteristics  needed  if  non-constant  cross- 
section  rods  are  to  be  studied. 

The  paper  ooncludes  with  the  proof  of  three  theorems  which  give  some 
insight  into  the  behaviour  of  the  solution  of  the  generalized  Riccati  equation 
in  terms  of  simple  conditions  placed  on  the  coefficients  and  on  the  initial 
condition  for  the  acceleration  wave.  It  is  shown  that,  in  general,  the 
intensity  of  an  acceleration  wave  will  become  infinite,  leading  to  shock  wave 
formation,  after  propagating  only  a  finite  distance  along  the  rod.  An 
interesting  new  result  which  arises  out  of  this  analysis  is  the  existence  of  a  j 

threshold  for  the  initial  wave  intensity.  This  has  the  property  that  the 
closer  the  initial  intensity  is  to  this  threshold,  the  further  will  the  wave 
propagate  before  a  shock  wave  forms,  provided  only  that  the  initial  intensity 
exceeds  this  positive  threshold. 
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ACCELERATION  WAVE  PROPAGATION  IN  HYPRELASTIC  RODS 
OF  VARIABLE  CROSS-SECTION 

A.  Jeffreyt 

1.  introduction 

The  nonlinear  elastic  material  studied  in  this  paper  is  one  in  which  the 
stress  potential  £  is  dependent  only  on  the  deformation,  or  displacement, 
gradient.  It  is  thus  a  material  which  belongs  to  the  class  known  as 
hyperelastic  [1],  and  we  shall  consider  a  special  problem  concerning  one¬ 
dimensional  acceleration  wave  propagation,  so  that  E  will  depend  only  on  p, 
the  displacement  gradient  in  the  direction  of  propagation. 

Many  authors  have  considered  both  static  and  dynamic  problems  for  such 
materials,  of  whom  we  mention  only  Antman  [2]  and  Antman  and  Jordan  [3]  who 
studied  the  Kirchhoff  problem  for  nonlinearly  elastic  rods  and  qualitative 
properties  in  general,  Jeffrey  and  Teymur  [4]  and  Jeffrey  and  Suhubi  [5]  who 
considered  shock  wave  formation  and  acceleration  wave  propagation  through 
periodically  layered  media,  and  Antman  and  Liu  [6]  who  made  a  detailed  study 
of  travelling  waves  in  hyperelastic  rods  which  were  permitted  various  forms  of 
deformation.  In  what  follows  we  show  how,  when  a  variable  area  of  cross- 
section  occurs  in  a  rod,  the  effect  of  the  geometry  is  to  produce  a  transport 
equation  for  the  acceleration  wave  intensity  which  is  more  general  than  usual. 

In  Sections  2  and  3  of  this  paper  our  purpose  will  be  to  show  that  the 
transport  equation  for  acceleration  wave  propagation  along  the  axis  of 
symmetry  of  a  hyperelastic  rod  of  slowly  varying  cross-section  is  a 
generalized  Riccati  equation,  and  to  deduce  the  form  of  its  coefficients. 

This  is  a  generalization  of  the  transport  equation  studied  by  Bailey  and  Chen 
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[7,8,91  which  was  a  variable  coefficient  Bernoulli  equation  known  to  apply  in 
diverse  situations. 

Since  the  Bernoulli  equation  may  be  solved  exactly,  Bailey  and  Chen  were 
able  to  make  precise  estimates  concerning  the  growth  and  decay  properties  of 
solutions  in  terms  of  the  coefficients  and  the  initial  condition.  Such 
precise  estimates  are  not  possible  in  this  case,  so  in  Section  4  three 
representative  theorems  are  proved  concerning  the  unbounded  nature  (escape)  of 
the  acceleration  wave  intensity  s  after  propagation  along  only  a  finite 
length  Xn  of  such  a  rod.  It  is  shown,  for  example,  that  some  hyperelastic 
materials  and  cross-sectional  area  variations  always  lead  to  the  escape  of  s, 
whereas  in  others  a  threshold  a*  exists  for  s  with  the  property  that  the 
escape  can  only  be  certain  to  occur  if  s  >  0*. 


2.  Formulation  of  Problem 

We  shall  consider  a  homogeneous  hyperelastic  medium  [1]  in  the  form  of  a 
semi-infinite  rod  of  variable  cross-section,  which  has  the  stress  potential 
E.  Furthermore,  we  shall  assume  that  the  cross-section  is  everywhere 
symmetrical  relative  to  the  axis  of  the  rod  and  to  two  fixed  mutually 
orthogonal  axes  that  are  normal  to  it,  but  that  the  cross-sectional  area  is  a 
slowly  varying  function  of  distance  measured  along  the  rod.  The  material  will 
be  assumed  to  be  incompressible. 

Let  us  now  relate  the  Euler ian  coordinates  x^  (k  =  1,2,3)  to  the 
Lagrangian  coordinates  X^(k=  1,2,3)  by 


x  =  u(X,t)  +  X,  x1  =  x,  X1  =*  X, 


x,  =  X. , 

3  3 


(2.1) 


where  u(X,t)  is  the  displacement  of  a  particle  along  the  rod  relative  to  the 
natural  axis  along  which  X  is  also  measured  [10].  The  orientation  of  the 
X2  and  X^  axes,  which  are  mutually  orthogonal  and  lie  in  a  plane  normal  to 
the  X  axis,  is  taken  so  that  the  cross-section  is  symmetrically  disposed 
relative  to  them. 


In  terms  of  these  coordinates,  the  displacement  gradient 

&u 

p  "  "Sx  "  ux  ' 

while  the  velocity  of  a  particle  displaced  by  an  amount  u  is 


The  stress-tensor  is  given  by 


3u 

*t  "  Ut  * 


T11  *  T  "  3p  ' 


(2.2) 


(2.3) 


(2.4) 


where  because  of  the  hyperelasticity  E  ■  E(p). 

Since  the  particle  displacement  u  will  be  taken  to  be  small,  we  shall 
regard  the  cross-section  area  variation  solely  as  a  function  of  X,  and  denote 


t 

m, 
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it  by  S(X).  We  suppose  the  rod  to  be  stressed  and  left  in  equilibrium  prior 
to  the  propagation  of  a  disturbance.  Then  continuity  of  stress  normal  to  the 
(X2,X3)-plane,  coupled  with  the  area  variation  iteself ,  shows  that  the  initial 


displacement  p^(X)  will  be  given  by  solving 


;  .  siil  r«)  t 

S(0>  * 


(2.5) 


*  S(0 )  v3p'  * 

where  T  »  T(0)  is  the  initial  constant  stress  at  the  plane  boundary  X  -  0, 
which  we  take  to  be  the  end  of  the  rod  at  which  the  cross-sectional  area  is 


S(0>. 


The  equation  of  motion  is  simply 

3(TS) 

3X  " 


SP  It  ' 


(2.6) 


where  P  is  the  constant  density.  Then,  by  virtue  of  equation  (2.4), 

3t  3t  3^  32E  3p 

5x  "  3p  3x  "  3p2  3x* 

If  we  now  introduce  the  quantity  c,  which  has  the  dimensions  of  a 
velocity,  through  the  definition 


r\  d2Iy6  i2Z 

c  -  c(p)  -  — -J  ,  with  — -  >  0  , 

P  3p2  3p2 


(2.7) 


the  equation  of  motion  (2.6)  becomes 


3v  2  3p  T  d  . .  . 

^  -  c  -  p  a  {tn  s)  "  0  ' 


(2.8) 


where  T  -  T(X)  is  the  stress  as  a  function  of  X,  which  is  influenced  by 
both  the  propagation  of  the  wave  and  by  the  change  of  equilibrium  stress 
caused  by  the  variation  of  the  cross-sectional  area. 


Equality  of  mixed  derivatives  then  leads  to  the  additional  equation 

3p  3v 
3t  "  3X  * 


(2.9) 


Thus  the  one-dimensional  theory  for  the  propagation  of  disturbances  along  the 
rod  is  governed  by  the  equations  (2.8)  and  (2.9).  The  last  term  in  (2.8) 
represents  the  effect  of  the  cross-sectional  area  change,  and  when  S  -  const. 


-4- 


these  equations  reduce  to  those  considered  elsewhere  in  a  somewhat  different 
context  [4,5].  This  is  because  the  one- dimensional  model  does  not  distinguish 
between  wave  propagation  normal  to  the  bounding  plane  of  a  half-space  and  wave 
propagation  along  a  rod  of  constant  cross-section. 
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3.  The  Transport  Equation  For  The  Acceleration  Wave  Amplitude 

We  shall  suppose  that  a  wave  propagating  along  the  the  rod  in  the 
positive  X  direction  starts  from  the  end  of  the  rod  corresponding  to  X  * 

0  at  time  t  =*  0,  and  that  it  advances  into  the  equilibrium  state  determined 
by  (2.5).  We  thus  assume  that: 

(i)  p  and  v  are  continuous,  with  v(X,t)  ■  0  ahead  of  the  wave, 

(ii)  the  first  and  second  derivatives  of  p  and  v  experience  at  most  a 
jump  discontinuity,  so  that  the  wavefront  which  propagates  is  an 
acceleration  wave. 

Let  us  now  generalize  an  argument  used  by  Gurtin  [11]  in  connection  with 
water  waves,  and  show  that  it  may  also  be  used  in  this  case,  though  we  shall 
see  that  the  transport  equation  which  results  differs  significantly. 

In  what  follows  we  denote  by  a  superscript  minus  sign  the  value  of  a 
quantity  immediately  behind  the  advancing  wavefront,  that  is  to  say  in  the 
disturbed  region,  and  by  a  superscipt  plus  sign  the  corresponding  value  of 
the  wavefront  in  the  equilibrium  state. 


We  then  conclude  from  (i)  above  that 


p+  and  v  =  v+  =  0  , 


(3.1) 


while  from  (2.8)  it  may  be  seen  that  the  time- independent  solution  p+  must 


satisfy 


co  If +  r  h  <*n  s>  - 0 ' 


(3.2) 


where  the  suffix  zero  signifies  the  form  of  the  function  on  the  wavefront. 


Differentiating  v  parallel  to  and  immediately  behind  the  wavefront  we 


obtain 


3v  dt  3v 
3X  dX  3t 


(3.3) 


to  which  result  may  be  supplemented  the  form  taken  by  (2.9)  immediately  behind 
the  wavefront 


(3.4) 


3p  3v 

3t  mJT  * 

Next,  differentiating  p  »  p+  in  similiar  fashion,  and  noting  that  p+  is  a 


function  only  of  X,  we  find 


£e_  +  3g^ 

3x  dX  3t  3x  ' 


which  by  virtue  of  (3.2)  becomes 

3p  dt  3p  T0  d  . 

st  +  ffl  if"  +  ~1  3c  (£n  S)  "  0  ‘ 

Pc, 


(3.5) 


(3.6) 


dX  -)2  2 

—  ~  c0  * 


(dX) 

Mt; 


This  must  also  be  supplemented  by  the  form  taken  by  (2.8)  immediately  behind 
the  wavefront 

I?" '  co  if  ■  p3,  h  (Zn  s)  *  °*  (3*7) 

We  are  now  in  a  position  to  interpret  dX/dt,  the  speed  of  wavefront 
propagation,  in  terms  of  the  material  constants.  Combining  equations  (3.3), 
(3.4),  (3.6)  and  (3.7)  and  solving  for  dX/dt  we  find  that  in  fact 

(3.8) 

Thus  the  quantity  c  introduced  in  (2.7)  is  merely  the  propagation  speed  of 
the  disturbance  (acceleration  wave)  characterized  by  (i)  and  (ii)  above.  We 
notice  that  TQ  =  TQ(X)  and  cQ  =  cQ(X),  since  the  initial  displacement  given 
by  solving  (2.5)  is  pQ  =  pQ(X). 

Hereafter  we  confine  attention  to  the  situation  immediately  behind  the 
wavefront.  Differentiating  equations  (2.8)  and  (2.9)  partially  with  respect 
to  t  and  X,  respectively,  and  eliminating  3  p  /3t3x  gives 


32v~ 

3t2 


-  c. 


2  32  v~ 
ax2 


'  J  fc  ««■•>-  • 

3p 


1  3 


(3.9) 


Now  it  follows  from  (3.8)  that  on  the  wavefront  3/3t  =  cQ  d/dX,  which 
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leads  at  once  to  the  Identity 

aV  *JV 
at2  *  C°  ® 

It  is  also  true  that  on  the  wavefront 

aV  .  *"1' 

.2  "  dX  ' 


3x‘ 

so  combining  (3.10)  and  (3.11)  brings  us  to  the  final  identity 


a2  - 
3  v 


„  *2  - 
2  3  v 


3t2  ’  C°  3X2 


d(vfc) 


2  d(V 


'0  dX  C0  dX 


(3.10) 


(3.11) 


(3.12) 


Setting  3p  /3x  *  s,  the  acceleration  wave  intensity,  and  combining  (3.9) 
and  (3.12),  brings  us  to  the  result 


ds  .2 

—  »  -Ps  +  Bs  +  Y  , 


where 


P(X) 


B(X)  = 


Y(X) 


_1 _ 

r  3c 

_ 2.  f. 

2pco 

L  p  ^ 
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rli£i 

4Pco 

a  30 
3p 

V 


3p 


S(X)) 


(3.13) 


(3.14) 


(3.15) 


2Pc: 


2  2  T 

2c  — •  (An  S(X))  +  (^-|)  (An  S(X))}2  .  (3.16) 

.  dX  0  3p  0  J 


Equation  (3.13)  is  the  transport  equation  for  the  acceleration  wave 
intensity,  and  it  has  been  expressed  in  this  form  as  this  is  an  obvious 
generalization  of  the  canonical  form  introduced  by  Bailey  and  Chen  [7,8,9]. 
However,  on  account  of  the  geometrical  effect  introduced  by  the  change  of 
cross-sectional  area,  (3.13)  is  now  a  generalized  Riccati  equation  [12], 
rather  than  the  Bernoulli  equation  found  by  Bailey  and  Chen. 

We  conclude  this  section  with  some  general  remarks  about  the  coefficients 
P,  B,  and  Y.  As  observed  by  Bailey  and  Chen,  the  coefficient  B  is  soley 
determined  by  the  material  behaviour,  but  P  and  the  new  coefficient  Y 
involve  both  the  material  behaviour  and  the  geometry.  Now  >  0,  and  it  has 


already  been  observed  in  (2.7)  that  for  c  to  be  real  it  is  necessary  to  have 


32£/3p2  >  o.  However  the  sign  of  33T/3p3  will  be  determined  by  the 

particular  material  that  is  involved. 

Some  special  cases  of  area  variation  arise  which  are  worthy  of  note. 
d2 

If  S  is  such  that  — -  (in  S)  >0  the  graph  of  In  S  is  concave-up,  and  it 

<K2  d 

will  be  increasing  if  —  (in  S)  >  0  and  decreasing  if  —  (in  S)  <  0  . 

dX  „  dx 

cl 

Conversely,  if  S  is  such  that  — —  ( An  S)  <  0  the  graph  of  £n  S  is 

dX2  d 

con cave- down ,  and  it  will  be  increasing  if  —r  ( in  S)  >0  and  decreasing  if 

aX 

(*n  8)  <  0. 

When  S  =  const,  both  y  and  Y  vanish  and  8  =  =  const.  Equation 


(3.13)  may  then  be  integrated  to  give 


s(X) 


i(0) 


1-0os(O)X 


(3.17) 


where  s(0)  is  the  initial  acceleration  wave  intensity  when  X  *  0.  In  the 


event  that  B0s(0)  >  0  the  acceleration  wave  intensity  will  become  infinite 
after  it  has  propagated  a  finite  distance  XB  ,  which  we  shall  call  the  escape 
distance.  In  this  case  this  is  given  by  [see  4) 


x®  =  V8  s(0)  .  (3.18) 

When  $Qs(0)  <  0  the  acceleration  wave  intensity  will  simply  decay  to  zero  as 
the  wave  propagates.  The  position  X  »  XB  may  be  interpreted  as  the  point  at 
which  a  shock  first  forms,  but  we  shall  not  pursue  this  aspect  further  and 
refer  instead  to  Chen  [9] . 

The  transport  equation  (3.13)  can  only  be  used  to  study  general  wave 
2  2  3  3 

propagation  if  3  E/3p  and  3  E/3p  are  known,  and  the  simple  case  of  a 
constant  cross-section  rod  of  fixed  length  L  may  be  used  to  provide  this 
information.  If  the  sound  speed  c(p)  in  such  a  rod  is  measured  for 
different  displacement  gradients  p,  then  3  £/3p  follows  by  use  of  (2.7). 
Similiarly,  if  for  such  a  rod  of  length  L  the  initial  intensity  s(0)  and 

3  3 

transmitted  intensity  s(L)  are  measured  for  different  p,  then  3  I/3p 


■9 


follows  by  solving  (3.17)  for  to  get 


B. 


s(L)  -  s(0) 
Ls(  0 ) 


(3.19) 


and  then  using  (3.15).  This  pre-supposes ,  of  course,  that  the  length  L  is 
sufficiently  short  to  avoid  escape  of  the  intensity  s  before  the  end  of  the 


rod  is  reached. 


Escape  of  Solutions  to  the  Transport  Equation 
In  this  final  section  we  offer  some  results  concerning  the  conditions 
necessary  for  the  escape  of  the  solution  of  the  transport  equation  (3.13). 
These  are  obtained  by  making  simple  assumptions  about  the  coefficients 
a,  B  and  Y  of  equation  (3.13),  and  then  using  elementary  comparison 
equation  methods.  No  systematic  study  of  the  qualitative  properties  of  the 
transport  equation  has  been  attempted,  since  our  objective  is  merely  to 
establish  that,  in  general,  an  acceleration  wave  in  a  hyperelastic  material 
will  have  a  finite  escape  distance.  However,  in  addition,  we  also  establish 
the  existence  of  thresholds  o*  >  0  and  o^  <  0.  These  have  the  property 
that  in  certain  circumstances  escape  only  occurs  at  an  infinite  distance  when 
the  initial  acceleration  wave  intensity  s(0)  is  such  that  s(0)  +  , 

while  when  s(0)  t  a~  the  solution  s  is  monotonic  increasing  with 
lim  s(X)  >  a~  . 

X-M» 

We  now  prove  three  theorems,  the  first  two  being  related,  and  begin  by 

supposing  Y(X)  >  0  for  X  >  0,  for  then 

2  2 
-Us  +  Bs  +  Y  >  -Us  +  Bs  . 


Now  consider  the  equation 


ds 

dX 


Us  +  Bs 


2 


+  Y  , 


(4.1) 


and  the  comparison  equation 

~  -  go  +  Bo2  ,  (4.2) 

aX 

with  s(0)  =  o(0)  >  0,  where  we  now  assume  U(x)  <  0  ,  B(X)  >  0  for  X  *  0. 
Then 


d_ 

dX 


(s  -  O) 


and  it  follows  directly  that 


*  Y(0)  >  0  , 

X  =  0 


d_ 

dX 


(s  -  0)  >  0 


for  X  >  0  , 


unless  the  derivative  becomes  unbounded  for  some  finite  escape  distance 


(4.3) 

+ 

X  . 

CO 


-11- 


■tr 


i 


This  shows  that  s  >  0  as  long  as  the  derivative  in  (4.3)  is  defined.  If  o 

"“+  + 

has  a  finite  escape  distance  XM,  then  so  also  has  s,  say  X,,,  ,  where 

+  "f  ■>»f  +  •»«+ 

0  <  X^  <  Xw  showing  that  X^  provides  an  upper  bound  for  X,,.  To  find  X,, 
we  may  use  the  results  of  Bailey  and  Chen  [7,  Theorem  4.3(i)  or  9,  Theorem 
3.2.4(i)]  in  connection  with  the  comparison  equation  (4.2)  with 
0(0)  -  s(0)  >  0  . 


The  result  asserts  that  if  s(0)  >  a  ,  where 

“  X 

a  =  1/C/  8(X)  exp(-  /  p(T)dx)dx}  , 

0  0 

then  there  exists  a  finite  escape  distance  XQ  >  0  such  that 


(4.4) 


/  8(X)  exp(-  /  U(T)dT)dX  -  1/s(0). 
0  0 


(4.5) 


We  conclude  that,  when  a  finite  X^  exists,  there  is  some  X,,,  with 
0  <  X1  <  *1  such  that 

OB  OB 

lim  s(X)  »  +°°  . 
x*x* 

OB 

The  second  part  of  the  theorem  by  Bailey  and  Chen  cannot  be  used  to 
establish  the  decay  of  s(X)  to  zero  as  X++*  when  0<s(0)<a.  This  is 
because  we  are  working  with  the  inequality  s  >  0  and  the  decay  of  0  does 
not  imply  the  decay  of  s. 

Let  us  now  consider  equation  (4.1)  and  the  comparison  equation  (4.2)  when 
U(X)  <  0,  0(X)  <  0,  Y(X)  <  0  and  s(0)  -  0(0)  <  0  .  Then  the  same  form  of 
reasoning  as  before  again  establishes  that  escape  is  possible  for  o  at  some 
escape  distance  XM  >  0  ,  determined  by  the  results  of  Bailey  and  Chen  when  a 

threshold  condition  on  s(0)  is  satisfied,  and  that  this  time 

lim_0(X)  -  . 

X+x“ 

00 

Since  the  stated  conditions  ensure  that  s  <  0  ,  we  conclude  that  when  0  has 

an  escape  distance  X,,,  ,  then  so  has  s,  say  X_  ,  where  0  <  X_  <  X_  ,  and 
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that 


lint  s(X)  -  -® 

x+x" 

w 

For  the  same  reason  as  before,  no  conclusion  may  be  reached  from  the 
second  part  of  the  theorem  by  Bailey  and  Chen  concerning  the  decay  of  s  to 
zero  when  the  threshold  condition  on  s(0)  is  not  satisfied  and 
-a  <  s(0)  <  0  . 

We  have  proved  the  following  result. 

Theorem  1 

Consider  the  differential  equation 

~  -  Ps  +  Bs2  +  Y  , 

aX 

where  U,  B  and  Y  are  integrable  functions  of  X  on  every  finite  subinterval 
of  [0,®)  ;  and  let 

®  X 

a  a  1/  {  /  | B(X) |  exp(-  /  U(T)dT)dX>  . 

0  0 


fp<X)  <  0,  3(X)  >  0,  Y(X) 
<  0,  3<x)  <  0,  Y(X) 


with  |s(0) I  >  a  , 


mu  if  /  '■>“(  "  “  I  >->"l  '  “I  •»■<»/  >  ol  .  fs(0)  >  o\ 

Then/  if  \ n/V)  ✓  a  o/v)  >>  a  v/v\  ^  qj  |s(0)  <  0J 

>  o' 

there  exi§£s  a  unique  finite  escape  distance/  )  for  s(X),  with  upper 

O  > 0 

bounds  ^  y,  such  that 

<  , 

X* 

®  X 

/  B(X)  exp (-/  M(T)dT)dX  =*  1/s(0)  , 

0  0 


and 


flim  s(X) 

J 

|  lim  s(X) 

Lxfx- 


+® 


Now  suppose  constants  p^,  6^  exist  such  that  P(X)  <  p^  <  0, 
8(X)  >  BQ  >  0,  and  suppose  also  that  Y(X)  >  0  and  s(0)  >  0.  Then 
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-  Us  +  3s2  +  Y  >  -  WQs  +  Bq82  , 

so  that  considering  equation  (4.1)  and  the  constant  coefficient  comparison 
equation 

do  a  2 

«"'V  +  V  (4*6> 

with  0(0)  -  s(0)  >  0  we  see,  as  before,  that  s  >  0.  Applying  Theorem  1  to 
(4.6)  we  find  <*  ■  0,  so  that  in  this  case  an  escape  distance  X*  always 

exists  for  o,  and  hence  for  s.  A  trivial  integration  of  (4.6),  or  of  the 

's*+ 

condition  determining  x^  in  Theorem  1,  shows 

X»  =  An^s(0)-Xo^  >  0  '  (4*7) 

where  X^  =  U^/B^  >  0.  Hence,  if  X*  is  the  escape  distance  for  s,  we  have 

0  <  X*  <  X*  ,  where 

lim  s(X)  -  +*  . 

x+xl 


A  similar  form  of  argument  applies  when  W(X)  <  jjq  <  0,  3(X)  <  Bq  <  0, 
Y(X)  <  0  and  s(0)  <  0,  only  now  the  solution  s  escapes  negatively  at  some 
escape  distance  X~  where  0  <  X~  <  x~  and  x£  is  still  given  by  equation 
(4.7).  We  have  thus  arrived  at  our  next  result. 

Theorem  2 

Let  s  satisfy  the  differential  equation 

ds  a  2 

~  -  Us  +  Bs  +  Y  . 


Then  if  constants  U^,  B^  exist  such  that 

fyU)  <  uQ  <  o,  B(x)  >  B0  >  (A  JY(X)  >  o,  s(0)  >  o^l 

\uu)  <  UQ  <  0,  S(X)  <  B0  <  oj and  \y(X)  <  0,  s(0 )  <  0  J 


\_Y  (X) 

follows  that  there  exists  a  finite  escape  distance 


for  X  >  0,  it 


lim  s(X) 
X+X* 

OB 

|  lim  s(X) 

x+x" 


I  x*  >  o 


X.  >  o 


for  s  such  that 


_«•/'  w^ere  an  upper  bound 


on 


is  provided  by 
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with 


0 


0 


X 


0 


VBo  <  °\ 

v»o  *  °j 


For  our  final  theorem  we  begin  by  supposing  constants  a^,  3^,  Yfl  exist 
such  that  M(X)  <  <  0,  3<X)  <  Bfl  <  0  and  YQ  <  Y(X)  <  0  for  X  >  0,  and 

that  s  satisfies  equation  (4.1)  with  an  appropriate  initial  condition  s(0) 
which  ensures 

-Us  +  Bs2  +  Y  >  -W0s  +  0Qs2  +  YQ  >  0  .  (4.8) 

Clearly  condition  (4.8)  will  be  true  at  X  *  0  if  |s(0)|  is 
sufficiently  large,  thereby  ensuring  that  when  comparing  equation  (4.1)  with 
the  comparison  equation 

s  -  -V  *  “o'”2  +  \  ■  ,4-9> 

the  inequality  will  remain  true  for  X  >  0  so  that  s  >  o  .  Hence  the  escape 
of  o  will  ensure  the  escape  of  s. 

Now  the  right  hand  side  of  (4.9)  will  be  positive  provided  o  does  not 
lie  between  its  two  zeros 

•5  *  <*0  *  <“2.  -  4SoVV2,/26o  • 

Consequently,  inequality  (4.8)  will  be  true  for  X  >  0  ,  so  that  s  >  o, 
provided  either  (i)  s(0)  »  0(0)  >  0*  >  Jb,  or  (ii)  s(0)  *  0(0)  <  0fl  <  0. 
Suppose  (i)  is  true  and  write  (4.9)  in  the  form 


JL  31  .  v2 

s0  a 


2 

a  , 


(4.10) 


where  Y  •  o  -  (PO/20Q)  and  a2  -  (V  /20o)2  "  **0^0*  or»  e<Iuival*ntly» 


i 

( 

i 


a  *  o*  -  (W0/28q)  with  a  >  0.  Then  (4.10)  may  be  integrated  to  give 

-  *■  ■ 

where  YQ  =*  0(0)  -  (yo/28Q). 


(4.11) 


This  may  be  solved  for  0  when  we  find 

U 

a  ■ 


.  1  +  A  exp(23.aX) 

o  r _  Q  ) 

2$Q  *M  -  A  exp(20oaX)J  ' 


where 


Va  a(0)  -  a0 


V 


and  0  <  A  <  1. 


0(0)  -  0„ 


(4.12) 


Inspection  of  (4.12)  shows  that  escape  of  o  occurs  at  X,,,,  when 
1  =  A  exp(2PQa  X,,,).  This  leads  to  the  required  upper  bound  x,,,  for  the 


escape  distance  for  s 


230a 


o(0)  -  an 

in( - M  >  0  . 

0(0)  -  oft 


(4.13) 


The  significance  of  the  threshold  o*  is  seen  by  examining  result 
(4.13).  The  closer  the  initial  value  s(0)  =  o{0)  is  to  o*,  the  further  the 
comparison  solution  propagates  before  the  escape  of  o  occurs. 


Repetition  of  the  same  form  of  argument  for  case  (ii)  leads  to  the 


results 

HQ  r1  +  B  exp(2BQbX)^ 

0  *  2iT  +  b  1  -  B  exp(20QbX)  ' 
where  b  *  (11^/2  BQ)  ”  °0  w*th  ^  >  0,  and 


Y0  -  b  0(0)  -  0+ 

B  *  — - —  *  1  . . .  with  B>1.  (4.15) 

Yo  +  b  o(0)  -  o“ 


Inspection  of  (4.14)  shows  that  since  B  >  1  and  S^b  >  0,  there  can  be 
no  escape  of  0  as  X  ♦  +».  In  fact,  as  X  ♦  +•  so  a  *  (Wo/20o)  -  b  •  o“, 
showing  that  as  X  ♦  +“  s  must  exceed  the  value  o^. 

We  have  thus  proved  our  final  result. 
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Theorem  3 


Let  s  satisfy  the  differential  equation 

H  «  -us  +  0S2  +  Y  , 

subject  to  the  initial  condition  s  =  o( 0),  where  U(X)  <  <  0, 

0(X)  >  <  Y(X)  <0  for  X  >  0,  and  set 

o*  *  [u  +  (u^  -  4B  Y  )  ^2  i/2g  . 

0  0  1  0  00  J/  0 

Then  it  follows  that: 

(i)  if  0(0)  >  0+  the  solution  s  will  escape  at  a  finite  escape 

+  *vf  + 

distance  X^  >  0,  where  an  upper  bound  Xw  on  is  provided  by 


C260<  -  “o' 


£n(- 


0(0)  -  or 


0(0)  -  O 


f) 


and 


lim  s(X)  -  +®  ; 

x+x. 


(ii)  if  0(0)  <  0q  the  solution  s(X)  is  monotonic  increasing  and 

lim  s(X)  >  o“  . 

X-H* 

In  conclusion,  we  remark  that  other  similar  theorems  may  be  formulated, 
but  these  three  will  suffice  to  indicate  the  behaviour  of  the  solution  to  the 
transport  equation  (3.13)  in  some  typical  circumstances.  Here  again  we  will 
not  pursue  further  the  matter  of  shock  wave  formation  associated  with  the 
escape  of  the  acceleration  wave  intensity  s. 
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